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Abstract 

Polymer entanglements lead to complicated topological constraints 
and interactions between neighbouring chains in a dense solution or 
melt. Entanglements can be treated in a mean field approach, within 
the famous reptation model, since they effectively confine each indi- 
vidual chain in a tube-like geometry. In polymer networks, due to 
crosslinks preventing the reptation constraint release, entanglements 
acquire a different topological meaning and have a much stronger ef- 
fect on the resulting mechanical response. We apply the classical 
ideas of reptation dynamics to calculate the effective rubber-elastic 
free energy of an entangled rubbery network. We then compare the 
results with other theoretical approaches and establish a particularly 
close mapping with the hoop-model, with equally good description 
of experimental data. The present consistent reptation theory allows 
further development of dynamic theory of stress relaxation. 

PACS numbers: 

61.41.-|-e Polymers, elastomers, and plastics 
62.20.Dc Elasticity, elastic constants 
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1 Introduction 



Rubbery polymer networks are highly complex disordered systems. The sim- 
plest theoretical models consider them as being made of "phantom chains" , 
where each polymer is modelled by a three-dimensional random walk in space. 
To form a corresponding phantom network, the chains are crosslinked to each 
other at their end points, but do not interact otherwise, in particular, they 
are able to fluctuate freely between crosslinks. This has the unphysical con- 
sequence that the strands can pass through each other. If one tries to avoid 
this assumption, the theory is confronted with the intractable complexity of 
entanglements and their topological constraints. The mean field treatment 
of entangled polymer systems is a now classical reptation theory (de Gennes 
1979; Doi & Edwards 1986), which had a spectacular success in describing a 
large variety of different physical effects in melts and semi-dilute solutions. 
However, the parallel description of crosslinked rubbery networks has been 
much less successful. First of all, one has to appreciate a significant difference 
in the entanglement topology: in a polymer melt the confining chain has to 
be long enough to form a topological knot around a chosen polymer; even 
then the constraint is only dynamical and can be released by a reptation 
diffusion along the chain path. In a crosslinked network, any loop around a 
chosen strand becomes an entanglement, which could be mobile but cannot 
be released altogether. A number of other complexities arise from possible 
nematic interactions between rigid chain segments (Abramchuk et al. 1989; 
Bladon & Warner 1993) and from the coupling between imposed deformations 
and chain anisotropy known as the stress-optical effects (Jarry & Monnerie 
1979; Deloche & Samulski 1981, Doi et al. 1989). 

An early model of elastic response of entangled rubbers was developed 
by Edwards (1977): in tradition with the melt theory, it assumed that the 
presence of neighbouring strands in a dense network effectively confines a 
particular polymer strand to a tube, whose axis defines the primitive path, 
see Fig. [TJ. Within this tube, the polymer is free to explore all possible con- 
figurations, performing random excursions, parallel and perpendicular to the 
axis of the tube. One can show that on deformation the length of the prim- 
itive path increases. Since the arc length of the polymer is constant, the 
amount of chain available for perpendicular excursions is reduced, leading 
to a reduction in entropy and hence to an increase in free energy. However, 
the particular calculation in Edwards (1977) has a number of shortcomings; 
perhaps the main limitation is that one only looks at the entropy reduction 
associated with the overall change of primitive path contour length on de- 
formation, ignoring the essential mechanisms of local reptation and segment 
re-distribution between different tube segments. 
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Gaylord and Douglas (GD) (1990) developed a simple "localisation 
model" for rubbers, based on scaling arguments. Each strand segment is 
thought to be placed in a hard tube of a square cross-section. Assuming that 
the deformation affinely changes the dimensions of the tube, one can calculate 
the change in free energy of a network of chains, each confined in such a tube. 
Ball, Doi, Edwards and Warner (BDEW) (1981) chose a completely different 
approach by introducing the "slip-link model" : the effect of entanglements is 
not considered to be a permanent one, changing the environment of a single 
polymer strand as in the two models presented above, but rather leading to a 
local mobile confinement site, a link between two interwound strands, which 
is able to "slip" up and down along both strands. In a further development 
of this idea, Higgs and Ball (HB) (1989) adopted a similar approach, which 
is however mathematically much simpler: the entanglements localise certain 
short segments of a particular strand to a small volume. One can model this 
effect by describing a network strand as a free Gaussian random walk, which 
is, however, forced to pass through a certain number of hoops, which are 
fixed in space. We shall find that the results of HB are very close to ours, 
in spite of a number of significant differences in the physical model. This 
mapping gives confidence in the final expression for the rubber-elastic free 
energy and the role of network entanglements in it. 

In addition to the mentioned above, one can find a variety of other theo- 
retical models, some of which are discussed in the review article (Edwards & 
Vilgis 1988). For instance, the constrained junction fluctuation model (Flory 
& Erman 1982), which assumes that the entanglements primarily affect the 
fluctuation of the junction points. An overview of different models can be 
found for example in the reviews by Heinrich et al. (1995) and Han et al. 
(1999). All theoretical models, describing the macroscopic equilibrium elas- 
tic response of densely entangled rubbery networks, have a common purpose 
- to develop a physically consistent description depicting what one accepts 
as a correct coarse-grained molecular behaviour - but also accounting for a 
number of experimental results showing substantial deviations in stress-strain 
response from the ideal phantom-network result a = //(A — 1/A 2 ) (Higgs & 
Gaylord 1990). Up to now, no theory succeeded on both of these fronts. 

In our current work, we develop a consistent implementation of the clas- 
sical tube model to take into account the entanglement effects in crosslinked 
polymer networks within the same framework as in the polymer melt dy- 
namics. We particularly focus on the entropy of internal reptation motion 
and resulting re-distribution of chain segments along the tube in a deformed 
state. The reptation tube model of entangled network strands provides a 
more accurate description in the sense that it keeps track of the allocation of 
chain segment excursions in the tubes. In this way, our model closes an im- 
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portant gap among the existing models with some unexpected consequences: 
the model, although based on existing ideas, analyses rubber elasticity in a 
new way The calculations, however, reveal that our results are very sim- 
ilar to the ones of the HB model, which approaches rubber elasticity in a 
different way. After brief recollection of the principles of network theory on 
the example of the ideal phantom network in the following section, section 
3 introduces the model and its properties in some detail, and outlines the 
derivation of the full expression for rubber-elastic free energy. Section 4 ex- 
plores the properties of the full expression, examining different limiting cases, 
as well as its linear-response limit. We conclude by comparing the concepts 
and the results of this work with previous theories. 

2 Classical phantom chain network 

Before considering densely entangled rubber, we briefly review the well- 
known results of the phantom chain network theory, which provides the basics 
to most other theoretical models. 

Assuming that a single polymer performs a free random walk in three 
dimensions, one finds that the end-to-end distance Rq obeys a Gaussian 
distribution in the long chain limit. This result goes back far in history: 
one can review its derivation and consequences in the classical text on this 
subject (Doi & Edwards 1986). The distribution of Ro is given by 

where h is the monomer step length and N the number of steps of the chain 
trajectory. The entropic free energy of such a random walk, therefore, is given 
by the logarithm of the number of conformations with the fixed Rq and has 
the form f3F = — lnPo(-Ro) = nNW^o + const, where (3 = l/k B T the inverse 
Boltzmann temperature. At formation of the network, i.e. at crosslinking, 
each chain in the polymer melt obeys the distribution ([I]), which is then 
permanently frozen in the network topology. 

One then assumes that the network junction points deform affmely with 
respect to their initial positions Rq following the macroscopic deformation 
described by the tensor A; hence we can write R = XRq. Therefore, the 
deformation A alters the free energy of each strand. The change of entropic 
free energy per chain of the whole network can be calculated by the usual 
quenched averaging: 

(3F = -(\nP(R)) P{Ro) = ^Tr(A T A), (2) 
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where we have dropped an irrelevant constant, arising as the logarithm of 
the normalisation in ([I]). The overall elastic energy density, in the first ap- 
proximation, is simply (0) multiplied by the number of elastically active 
network strands in the system per unit volume, which is proportional to 
the crosslinking density: 

^eiast = Tr (A T A) , with \i = n ch k B T. (3) 

Three positive remarks have to be made in defence of this simple model of 
rubber elasticity. First, consider the crosslinking connecting the end points of 
different polymer strands. Since the positions of end points of a single chain 
fluctuate strongly, the junctions reduce these fluctuations and therefore alter 
the single chain statistics. However, in spite of an apparent complexity, this 
effect merely introduces a multiplicative factor of the form 1 — 2/0, where <ft is 
the junction point functionality (see for example the review by Flory (1976)). 
Secondly, one can assume that the deformation preserves the volume, since 
the bulk (compression) modulus is by a factor of at least 10 4 greater than 
the shear modulus, which is proportional to /i; this implies the constraint 
det A = 1- Thirdly, the quenched average in equation @ does not average 
over chains of different arc lengths, but the fact that the result is independent 
of arc length, generalises the result to apply for chains of arbitrary length, 
or even for a polydisperse ensemble of chains. 

As a result if one considers, for example, a uniaxial extension A (by incom- 
pressibility, the two perpendicular directions will experience an equal contrac- 
tion l/\/A), the elastic energy would take the form F e i ast = |yu(A 2 + 2/A). 
This then allows us to calculate the nominal stress, which is the force divided 
by the area of the undeformed cross section: 

9F f^ M f A \ 

Cphantom — ~qT - /i (A - — I (4J 

Commonly, one draws the stress-strain curves in the Mooney-Rivlin represen- 
tation, which plots the reduced stress function /* = a/(A — 1/A 2 ) against 1/A 
(see Fig. ^ in the discussion below). This representation therefore indicates 
the degree of deviation from the simple phantom chain network behaviour. 

3 Reptation theory of rubber elasticity 

Following the original ideas of Edwards (1977), we assume that each network 
strand is limited in its lateral fluctuations by the presence of neighbouring 
chains. Therefore each segment of a polymer only explores configurations in 
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a limited volume, which is much smaller than the space occupied by an ideal 
random coil. Hence, the whole strand fluctuates around a certain trajectory, 
a mean path, which is called the primitive path in the reptation theory 
(Edwards 1977). The most intuitive way to visualise such a trajectory is 
to imagine that the chain is made shorter and thus stretched between the 
fixed crosslinking points. The taut portions of the chain will form the broken 
line of straight segments between the points of entanglement, which restrict 
the further tightening. This primitive path can be considered as a random 
walk with an associated typical step length, which is much bigger than the 
polymer step length, as sketched in Fig. [TJ. The number of corresponding 
tube segments M is determined by the average number of entanglements per 
chain (the situation with no entanglements corresponds to M = 1). 

Effectively, the real polymer is confined by the neighbouring chains to 
exercise its thermal motion only within a tube around the primitive path. 
Note that all the chains are in constant thermal motion, altering the local 
constraints they impose on each other. Hence, the fixed tube is a gross simpli- 
fication of the real situation. However, one expects this to be an even better 
approximation in rubber than in a corresponding melt (where the success 
of reptation theory is undeniable), because the restriction on chain repta- 
tion diffusion in a crosslinked network eliminates the possibility of constraint 
release. 

To handle the tube constraint mathematically, we traditionally assume 
that the chain segments are subjected to a quadratic potential, restricting 
their motion transversely to the primitive path. Along one polymer strand 
consisting of N monomers of effective step length b, there are M tube seg- 
ments, each containing s m , m = 1, . . . , M monomer steps. We infer the 
obvious condition 

M 

Y, s m = N. (5) 

771=1 

In effect, one has two random walks: the topologically fixed primitive path 
and the polymer chain restricted to move around it - both having the same 
end-to-end vector Rq, between the connected crosslinking points. 

Each tube segment m can be described by the span vector A m , joining the 
equilibrium positions of the strand monomers at the two ends of each tube 
segments. The number of tube segments M (or, equivalently, the nodes of 
the primitive path) is a free parameter of the theory, ultimately determined 
by the length of each polymer strand and the entanglement density. 

Since the primitive path is a topologically frozen characteristic of each 
network strand, we shall assume that all primitive path spans A m deform 
affinely with the macroscopic strain: A' m = AA m . This is the central point 
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Figure 1: A polymer strand is surrounded by neighbouring chains, which 
effectively confine the strand to a tube. The tube segment m, with the span 
vector A m along its axis, contains s m monomer steps, where m runs from 1 
to M. 



in the model: the rubber elastic response will arise due to the change in the 
number of polymer configurations in a distorted primitive path. To evaluate 
the number of conformations, we look separately at chain excursions parallel 
and perpendicular to the tube axis, along each span A m . Effectively, this 
amounts to introducing a new coordinate system for each tube segment, with 
one preferred axis along A m . In this direction one recovers the behaviour of 
a free random walk, giving rise to a one-dimensional Gaussian statistics in 
the long chain limit. Note that only one third of the steps s m in this tube 
segment would be involved in such parallel (longitudinal) excursions. We 
therefore obtain for the number of parallel excursions in a tube segment m: 

To determine the number of perpendicular (transverse) excursions, one 
can introduce the Green's function for the perpendicular steps made by the 
chain. In effect, we consider a two-dimensional random walk, with a total 
number of steps (2s m /3), in a centrosymmetric quadratic potential. For each 
of these two perpendicular coordinates, the Green's function satisfies the 
following modified diffusion equation (see e.g. Doi & Edwards 1986): 

d b 2 d 2 q 2 \ 

ds ~ Jdxj + ~2 x2f ) G ^' x / ;s ) = 8{x f -Xi)6{s), (7) 

where the Xi and Xf are the initial and final coordinates of the random 
walk with respect to the tube axis and q determines the strength of the 
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confining potential. This is an essential parameter of many reptation theories, 
directly related to the tube diameter a. Many discussions can be found 
in the literature regarding, for instance, the possible dependence of a on 
entanglement density and the possible change in a on deformation. However, 
we shall find that the potential strength q does not enter the final rubber- 
elastic energy, in an approximation that we expect to hold for a majority of 
cases. Such a universality resembles the situation with phantom networks, 
where many chain parameters do not contribute to the final results. 

The equation (J7|) is very common in the physics of polymers and its 
exact solution is known. However, we only need to consider a particular 
limit qobs m 3> 1 of this solution, which is the case of dense entanglements 
(resulting in a strong confining potential) and/or of a large number s m of 
monomers confined in the tube segment. Outside this limit, that is, when 
the tube diameter is the same order as the arc length of the confined chain, 
the whole concept of chain entanglements becomes irrelevant. In the case of 
our present interest, that is, in the strongly confined limit, the solution has 
a particularly simple form, see (Edwards (1977), 



G m (xi,x f ] s m ) oc exp ( -7^( x i + x )) ~ g?oK 



Remembering that there are two coordinates describing the transverse ex- 
cursions, we obtain for the two-dimensional Green's function of the tube 
segment m: 

G m (r u 77; s m ) oc exp (-^q bs m ^ exp ("|^(^ 2 + , (9) 

where r% and 77 are the initial and final transverse two-dimensional coordi- 
nates. 

The total number of transverse excursions is proportional to the inte- 
grated Green's function: 

Wt oc / dn I dr f G 



f^mK'ii 'ft s m)- 

Since the Green's function in the approximation (^) does not couple the 
initial or final coordinates to the number of segments s m , this integration 
will only produce a constant normalisation factor which can be discarded. 

Gathering the expressions for statistical weights of parallel and perpen- 
dicular excursions, we obtain the total number of configurations of a polymer 
segment consisting of s m monomers in a tube segment of span A m : 

W m = WlWt oc -^exp f-J?—A 2 m - \q bs m ) . (10) 
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Therefore, we find for the full number of configurations of the whole strand 
by summing over all M steps of the primitive path and integrating over 
the remaining degrees of freedom, the number of polymer segments confined 
between each node of the tube: 





W = / ds x ■ ■ ■ / ds M J] W m )S £ s m - N , (11 



where the constraint (^) on the polymer contour length is implemented by 
the delta-function. The statistical summation in (|TT|) takes into account the 
reptation motion of the polymer between its two crosslinked ends, by which 
the number of segments, s m , constrained within each tube segment can be 
changed and, thus, equilibrates for a given conformation of primitive path. 

Rewriting the delta-function as 8(x) = J dk e lkx , we proceed by finding 
the saddle points s* m which make the exponent of the statistical sum (|TT1) 
stationary. It can be verified that the normalisation factors 1 / a/s^ contribute 
only as a small correction to the saddle points 

s * m w {WQ^Tikj) ■ (12) 

The integral in ([□]) is consequently approximated by the steepest descent 
method. We repeat the same procedure for the integration of the auxiliary 
variable k, responsible for the conservation of the polymer arc length. The 



saddle point value k*, inserted back into (12) , gives the equilibrium number 
of polymer segments confined within a tube segment with the span vector 
A m : 

__ ^A m 

Z^ m =l 

where A m = |A m | is the length of the m-th step of the primitive path. 
We finally obtain the total number of configurations of one strand, confined 
within a tube whose primitive path is described by the set of vectors {A m }. 
The statistical weight W associated with this state is proportional to the 
probability distribution: 



W(A 1 , . . . , A M ) oc P({A m }) oc e-^ bN \ ' . (14) 



\ M-l 

V A/ A I 



This expression is a result corresponding to the ideal Gaussian P (Rq) in 
equation (0) for a non-entangled phantom chain. Note that the chain end-to- 
end distance Rq is also the end-to-end distance of the primitive path random 
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walk of variable- length steps: J2m=i A m = Ro- The probability distribution 
P({ A m }) for a chain confined in a reptation tube is reminiscent of a normal 
Gaussian, but is in fact significantly different in the form of the exponent 
and the denominator (note the modulus of the tube segment vector A m ). 

From equation fll4]) we obtain the formal expression for free energy of 
a chain confined to a tube with the primitive path conformation {A m }, 
which implicitly depends on the end-to-end vector Rq (the separation be- 
tween crosslink points): 

^ = 5^(f A.j'+^-DlngA.), (15) 



4 Free energy of deformations 

We now perform a procedure which is analogous to the one used to ob- 
tain equation (Q). In the polymer melt before crosslinking, the ensemble 
of chains obeys the distribution in (|H|) giving the free energy per strand 
flT5p. The process of crosslinking not only quenches the end points of each 
of the crosslinked strands, but also quenches the nodes of the primitive path 
A m , since the crosslinked chains cannot disentangle due to the fixed net- 
work topology. In our mean field approach, the tube segments described by 
A m are conserved, although they may be deformed by the strains applied to 
the network. For evaluating the quenched average, note that the statistical 
weight (|TJ]) treats the tube segments m in a symmetric way. This allows us 
to perform the explicit summation over the index m: 

(^(lr4 + (M - 1)ln (l>)) 

+(M-l)/lnf£ A m )V (16) 

where, in the second term, m and n ^ m are arbitrary indices of the primitive 
path steps. The brackets (• • •) refer to the average with respect to the weight 
P({A m }) given in ©. 

Furthermore, note that any affine deformation A transforms the vectors 
A m into A TO = AA m and, hence, their lengths A m = | A m | into A' m = |A_A m |, 
but leaves the quenched distribution P({A m }) unchanged. Bearing this in 
mind, we can evaluate the averages (|16D , leading to the free energy per chain 
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of the crosslinked network. The Appendix gives a more detailed account of 
how one evaluates the averages. The resulting average elastic energy den- 
sity also incorporates the density of crosslinked chains in the system by the 
constant \i = n ch k B T: 

F - = MfTT 1 ^ (17 > 

+ ^(M-l)^^(|Ai) 2 + MM-l)ln|A|, 
where the following notations are employed: 



|A| = — / dQ |A e| (11 



All J\e\=l 
1 



lnlAI = — dtt ln|A el. (19) 
An J\ e \=i ~ 

The notation TTT refers to the angular integration over a unit vector e. 

In the case of uniaxial deformation, where A takes a diagonal form with 
A" = A and X ± = l/y/X, the expressions (@|D and© can be calculated 
explicitly. The Appendix contains the formulae (equations ( |2"2| ) - (^)), which 
need to be inserted into fll7|) to give the full rubber-elastic energy density. In 
the small deformation limit, where A = 1 + e, expression ( |P7|) immediately 
leads to the Young's modulus E in F ext ~ \Ee 2 after noting that 

Tr(A T A) ~ 3 + 3 e 2 

TTT 2 „ 

A SS 1 + T^ 
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lnlAI w — e 2 . 
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For small simple shear, we can write the matrix A as \j — $ij + £UiVj, 
where u and v are two orthogonal unit vectors. As expected, the shear 
modulus G in F shear ps |Ge 2 is one third of the Young's modulus E: 

„ 1 /4 2M + 1 1. s 11M + 5N . s 

G=-E = u[ — + -(M-1) — . 20 

3 P V3 3M+1 5 V ; 3M + 17 v ; 

We point out that the expression ( fETf) is nearly identical to the rub- 
ber elastic free energy density obtained by Higgs and Ball (1989), although 
they start from a completely different set of physical and assumptions and 
mathematical framework. In HB, the entanglements make the strands to 
interact with each other only at certain points. This is modelled by forcing 
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the individual strands to pass through a certain number of hoops. They 
effectively divide the polymer strand into segments, each performing a three- 
dimensional phantom random walk. In this approach, the entanglements 
only exert constraints at points and allow the polymer to obey free Gaus- 
sian statistics between these interaction points. Nevertheless, the apparent 
identity of the final expression for the rubber-elastic free energy density is 
a comforting reflection of consistency in underlying physical concepts and 
the mathematical treatment of both models. The internal similarity of the 
models arises from the treatment of local chain reptation, equilibrating the 
distribution of its segments between the nodes of primitive path, before and 
after the deformation. 

From the expression (|T7D , we can recover the ordinary free energy of a 
phantom chain network by taking the case M — 1: -F e iast = |/iTr(A T A). 
This limit means physically that the polymer strand is placed in one single 
tube, tightly confined to the axis. Mathematically, a random walk in three 
dimensions with the total of N steps is equivalent to a random walk in 
one dimension along a given direction A, with N/3 steps, while the two 
perpendicular excursions in a tightly confining potential do not contribute to 
the global elastic response. This fact is the underlying reason why we recover 
the phantom chain network result by taking M — 1 in our model. Of course, 
in a parallel "hoop model" of HB the case M — 1 simply means that there 
are no constraints. 

On the other hand, as the number of tube segments M becomes very 
large, we obtain a rubber-elastic elastic energy of the form 



There are two ways to arrive at this limit of Af > 1 in a real physical 
situation: either the polymer melt is very dense, causing a high entanglement 
density, or the polymer chain is very long. In the latter case, the polymer 
strand experiences many confining entanglements along its path. 

Recall that the -F elast is the elastic energy density, which relates to the 
free energy per chain, such as the eq. (|15D, and is proportional to the density 
of elastic strands in the system: F elast cx \i = n cil k B T, where n ch is the 
number of crosslinked strands per unit volume. We can assume that in a 
melt, or a semidilute solution, the chain density is inversely proportional to 
the volume of an average chain, hence inversely proportional to the chain 
contour length: n c h oc 1/L. In the case of phantom network with a free 
energy F e i ast = |/uTr(A T A), one conclude that -F elast — > as the chains become 
infinitely long! This unphysical result reflects the fact that the phantom 
chain model assumes the entanglement interactions of the chains irrelevant. 




(21) 
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Clearly, this assumptions breaks down in the long chain limit, where we 
expect the entanglements to play a crucial role. 

This unphysical feature of phantom network model is overcome by our 
expression (f2T|). As the strands become longer, they will experience more 
entanglements, generating more primitive path nodes and confining tube 
segments. One can assume that the number of entanglements scales linearly 
with the strand length L: M oc L, in fact: M = N/N e with N e the char- 
acteristic entanglement length, essentially the average of s* m . Considering 
expression (|2lD , with n ch = 1/v monomcr iV (the inverse volume of the chain in 
a melt), we note that the corresponding elastic free energy F e i ast does not van- 
ish in the limit L — > oo. At M > 1 the constant plateau modulus becomes 
\xM = ksT / v monomer N e . As one expects, in an entangled polymer system 
there is no real difference between the network and the melt and they both 
have the same plateau modulus - with the corollary that the entanglement 
length N e must be much smaller in the permanently crosslinked network. 

5 Conclusion 

In this present work, we have analysed the behaviour of a polymer network 
in presence of entanglements, which are treated within a classical mean-field 
tube model. We found that this leads to an elastic energy which is identical to 
the one derived in the hoop model of HB (Higgs & Ball 1989). We claim that 
our model captures the physics of entanglements in a better way than HB, 
since the entanglements at a microscopic level do not so much localise the 
polymer at fixed points in space along its path, but rather impede the chain 
fluctuations on a length scale which is much larger than the size of monomers. 
This is due to the fact that the entanglements are caused by neighbouring 
strands, which likewise fluctuate. In this sense, our model provides a firmer 
ground of the theoretically known and experimentally tested results. 

Since our results for the isotropic rubber are identical to the ones by 
HB and theirs have been extensively compared with experimental data in 
a review by Higgs & Gaylord (1990), we do not need to include a similar 
comparison here. We only remark that the theoretically predicted curve for 
the reduced stress function /* agrees with the experimental data to a good 
level, see Fig. |[ 

Our model only describes the equilibrium response of a network to defor- 
mation. Shortly after applying the deformation, traditionally assumed affine 
on the level of primitive path segments, the network finds its chains far from 
the microscopic equilibrium. Each strand then redistributes its monomers 
between the deformed the tube segments, attributing more monomers to 
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Figure 2: The Mooney-Rivlin plot: the reduced stress function /* = 
er/(A — 1/A 2 ) plotted against 1/A. Data points from the analysis of ref. 
(Gaylord & Douglas 1990) are fitted by the result arising from full rubber- 
elastic energy fll7|) . Fitting indicates that M 3> 1 and the equation ( pi]) is a 
good approximation. 

some segments, less to others. This process proceeds via the local reptation, 
the sliding motion along the primitive path, and would be reflected by a 
time dependence of the variable s m . This number of monomers attributed 
to the tube segment m is changing, after the deformation, between a the 
initial value s* m in the Eq. (|13D and the corresponding saddle-point value at 
A' m = |A_A m |. By describing the associated relaxation process, our model can 
be naturally extended to describe the short-time viscoelastic response of an 
ideal crosslinked network (with no loops and dangling ends, leading to long 
relaxation of stress); the characteristic time of such a reptation dynamics is 
the Rouse time, since the dominant process is the monomer diffusion along 
the principal path. 

We appreciate many valuable discussions with S. F. Edwards. S. K. grate- 
fully acknowledges support from an Overseas Research Scholarship, from the 
Cambridge Overseas Trust and from Corpus Christi College. 

Evaluation of quenched averages (A m A n ) and 
(ln(EA m )) 

To evaluate the thermodynamic average (A m A n ), m,n = 1,...,N in the 
general case, one needs to average A m A n = |A m ||A n | with the probability 
distribution given by Eq. flT4]). For this purpose, one first has to find the 
normalisation M of the distribution, which can most easily be achieved by 
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introducing a new scalar variable u = Em=i to simplify the exponent of 
this distribution: 

m , exp (~Wn fem=i A r , 



(E A m ) M_1 



m=l 

M „ °? / M 

IJ J dA m J du5lu-Y<& 



— 3-« 2 

e 26 2 AT 



m=l n \ m=l / ^ 



m l „M-1 
u-Ai 



(4tt) m Jdu e -^JdA 1 Al JdA 2 A 



o 

-Ajf-2 



/ 



rfA M _ 1 A^„ 1 • (« - Ai - . . . - A M _x) 2 



In the last step, we introduced spherical coordinates for the variables A m , 
implemented the constraint u = Yl A m and used the fact that the variables 
A m are bound to be positive. The underlined expression is a function of 
u, which we call Im(u). Since the integrals only involve power functions, 
Im(u) itself is a power in u, whose order can be determined by counting the 
dimensions: 

Im{u) = -^u 3M ~\ 

After noting the recursive structure of Im(u), one can find a recursion relation 
for the coefficients X M : 

X M+1 = l-3M(3M + 1)(3M + 2)X M . 

Since Im=i(u) — u 2 , ah coefficients X M and hence all functions Im(u) are 
known. The remaining integration of u is a standard Gaussian integral. 

Having obtained the normalisation constant A/", we can proceed to calcu- 
late the averages (A m A n ). The calculations are very similar to the above one 
for Af; the only significant difference is the angular part of the integration 
which produces terms of the form §Tr(A_ T A_) and |AJ for the case m = n or 
m 7^ n respectively. 

For the logarithmic term, one in fact finds that only the angular integra- 
tion yields relevant terms: 
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1 f<X> e 262 at" 

du- 



M Jo 



u 



M-l 



u-Ai 



dA^f / dA 2 A^ 



— jw — js _ 

dA M -iA M _ 1 (it - Ai 



A M -i; 



^ II / dQ ™J M A i|Aei| + . . . + A M _i|Aej|f_i| 



+{u- Ax-...- A M -i)|Ae 



with the unit vector e m specifies the (arbitrary) orientation of the corre- 
sponding tube segment A m . Then one observes: 



In 



Ai|Aei| + . . . + A M _i|Ae M -i 
= ln[w] +ln[|Ae A/ |] +ln 



+ (u-A 1 



A 



M-l, 



lAejvfl 



M-l 



|Ae m | 



A r 



m=l \ 1= 



AeM I 



u 



small since «>|A r 



ln[|AeM|] + const. 



In case of uniaxial deformation, the evaluation of Tr(A T A), |A| and In |A|, 
cf. equations ( |T7D - (|1|) has been given elsewhere (Higgs & Gaylord 1990), 
but for completeness we state them here as well: 



Tr(A T A) 

|Af 
H\W) 



A2+ A 
A + 



1 
2 



In 



\3/2 _ V ^TT 



ln(A) - 1 



arctan \/A 3 — 1 



(22) 
(23) 

(24) 
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